Optical cloning of arbitrary images beyond the diffraction limits 
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Cloning of arbitrary images encoded onto the spatial profile of a laser beam onto that of a second 
beam is theoretically investigated. The two fields couple to an atomic lambda system in a coherent 
population trapping configuration. In particular, the case in which the probe and control fields 
are of comparable strength is considered. By considering more and more complex structures, we 
eventually find that our method is suitable to clone arbitrary images, which we demonstrated by a 
full numerical simulation of the propagation dynamics of both applied fields in the atomic medium, 
with the three letters "CPT" encoded on the initial control field profile. We find that the cloned 
structures have feature sizes reduced by about a factor of 2 compared to the initial images, consistent 
with a recent related experiment. 

PACS numbers: 42.50.Gy, 32.80.Qk, 42.65.-k 



I. INTRODUCTION 

Optical diffraction forms a fundamental limitation to 
the creation and detection of small images for conven- 
tional optical devices and microscopy. The reason is that 
any image with finite size can be considered as a group 
of different plane-wave components. Each component ac- 
quires a different phase shift during propagation either in 
free space or in medium. The resulting superposition of 
all wave components leads to diffraction, which means 
that the transmitted image will be distorted even after 
propagating only a few Rayleigh lengths |l| . From this in- 
terpretation it is clear that diffraction can be eliminated 
if conditions are tailored such that each plane- wave com- 
ponent acquires the same phase shift during propagation. 
To achieve this, it has for example been suggested that 
several specially-shaped beams with spatially-dependent 
phases such as Airy [2, Q: Bessel d, Q, and Mathieu 
beams p, lZ[, can propagate without diffraction in free 
space, since the phase differences between the different 
components are exactly compensated by the initially spa- 
tially dependent phases. 

It has also been recognized that optical diffraction 
can be greatly suppressed or even eliminated by using 
atomic coherence effects such as electromagnetically in- 
duced transparency |9l-ll5l| , coherent population trapping 
(CPT) [16l - |19| . or saturated absorption techniques [2y| in 
multilevel atomic systems. Most of these methods can be 
understood by noting that the spatially-dependent probe 
and coherent control fields couple to the atomic transition 
in such a way, that the probe field experiences a spatially 
dependent refractive index depending on the intensity of 
the control fields. As an example, a suitable spatially 
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dependent control field can optically induce a waveguide- 
like structure. When the probe field propagates in the in- 
duced waveguide inside the atomic medium, it is mostly 
confined in the waveguide due to total refiection, result- 
ing in the elimination of the diffraction for the probe. 

A more interesting and counterintuitive approach, 
which is fundamentally different from previous methods, 
has been developed by Firstenberg et al. [2l|, ^^ . They 
have theoretically proposed and experimentally verified 
that atomic motion and collisions can be utilized to elim- 
inate the paraxial diffraction of arbitrary image. The 
underlying physics is that the phase difference acquired 
during propagation for each plane- wave component of the 
image is exactly compensated by an additional phase 
shift induced by the atomic motion, thus leading to 
the elimination of diffraction. However, this method is 
plagued by strong absorption. 

Recently, Li et al. [23| experimentally demonstrated 
that the spatial shape of a control beam can be cast onto 
a weak probe beam via coherent population trapping in 
a three- level lambda atomic system. In their experiment, 
the transmitted intensity of the probe beam had a simi- 
lar spatial profile as that of the control beam, no matter 
what the input probe is. Moreover, the size of the trans- 
mitted probe beam was half of that of the diffraction- 
limited input probe. Manipulation of the susceptibility 
along the transverse direction is the key idea behind map- 
ping the spatial shape of the control beam onto the probe 
beam. It has been used for focusing, defocusing, and self- 
imaging SEHHi. 

It is important to note that the probe field is typi- 
cally treated as a weak field and the control field as a 
strong field, such that perturbation theory can be em- 
ployed to derive the linear effect of the atomic medium 
on the probe field propagation. There have been theo- 
retical and experimental studies where the probe is not 
necessarily weak |26l - l28l |. In this situation, the effect of 
the atomic coherence on the propagation dynamics of the 
control field need to be taken into account [231 ■ 

In this paper, we theoretically investigate the possibil- 



ity of cloning of an arbitrary image carried on an control 
field to another probe beam. The two control and probe 
beams are coupled to a three-level atomic lambda sys- 
tem to form a CPT configuration. We assume the two 
laser fields to be of comparable strength, such that per- 
turbation theory for the probe field is not valid any more 
to describe the effect of the atomic medium on the two 
fields. We start by calculating the susceptibilities includ- 
ing linear and nonlinear effects for both fields by solving 
the related density-matrix equations. As expected, we 
find that a spatially-dependent refractive index for the 
probe field is generated, structured by the spatial inten- 
sity profile of the control beam. In particular, the gener- 
ated structures enable one to transfer the transverse dis- 
tribution of the control field onto the transmission profile 
of the probe field. In order to study the full propagation 
dynamics, we then numerically solve the paraxial propa- 
gation equations for both fields by using a higher order 
split operator method [30|. We begin our analysis with 
a Gaussian control and a super- Gaussian probe field and 
observe the gradual mapping of the control field onto the 
probe field throughout the propagation. We in particular 
find that in the case of a strong probe field, the trans- 
mitted probe beam is focused more tightly by a factor of 
two compared to the weak probe field case. Next, we con- 
sider a control field with a spatial two-peaked Hermite- 
Gaussian profile, and demonstrated cloning of the profile 
onto the probe beam with feature size reduced by a factor 
of about 2.5. In order to verify that our method can serve 
as an universal tool for cloning of arbitrary image, we fi- 
nally simulate the three-dimensional light propagations 
for both fields, in which the spatial profile of the control 
field carries the three letters "CPT" . We show that also 
this structure can be cloned onto the probe beam which 
initially has a simple plane- wave profile, even though the 
control field is severely distorted throughout the propa- 
gation due to diffraction. Again, we observe a reduction 
of the feature size by a factor of about 2 in the probe 
field. 

The article is organized as follows. In Sec. II A, we 
introduce our theoretical model and then employ a semi- 
classic theory to describe the system. In Sec. II B, we dis- 
cuss the propagation equations for all fields in the strong 
field limit. In Sec. II C, we analytically calculate the 
atomic responses to the probe and control fields, includ- 
ing both linear and high-order nonlinear effects. In Sec. 

II D, we specify the spatially-dependent profiles for both 
fields that are used in our numerical simulation. In Sec. 

III A, we first show the differences in the spatial depen- 
dence of the absorption properties for the probe field in 
the weak and strong field limit, respectively. Then, we 
discuss the spatial refractive modulation which is depen- 
dent on the probe detuning. Finally, the refractive index 
for the control is investigated in both weak and strong 
field cases. In Sec. Ill B, we discuss our numerical results 
on the light propagation in several cases where different 
spatial profiles are chosen for both fields. Sec. IV pro- 
vides a brief summary of our paper. 
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FIG. 1. (Color online) Schematic setup for the cloning of 
arbitrary images encoded in the spatial profile of a control 
field onto the spatial profile of a probe beam. Both fields 
co-propagate through a rubidium vapor cell, and couple to 
the atoms on the |3) -f-^ |1) (probe field with frequency a;i) 
and 1 3) -f-*- |2) (control field with frequency 0)2) transitions, 
respectively. Our analysis includes the case in which both 
fields are strong. 



II. MODEL 

A. Equations of motion 

We consider a homogeneous cloud of ^^Rb atoms, with 
two laser fields coupled to the atoms such that a three- 
level A scheme as shown in Fig. [T] is of relevance. The 
probe field couples to transition |I) ^ |3), and the control 
field to transition |2) <-^ |3). The two co-propagating 
fields are defined as 

4(f,i)-e,f,(f)e-(-^*-'^^-)+c.c., (1) 

where £j (r) are the slowing varying envelope functions, ij 
the unit polarization vectors, tOi the laser field frequencies 
and kj is the wave numbers of fields. The index j € {1,2} 
denotes the probe or control field, respectively. 

The Hamiltonian of the system in electric dipole and 
rotating wave approximation is given by. 
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with Rabi frequencies of the probe and control fields de- 
fined as 
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The atomic transition frequencies and the corresponding 
dipole moment matrix elements are denoted by lo^j and 
dsj", respectively, and we label the radiative decay rate 
from state |3) to ground state \j) by 27^. The master 
equation for the density operator p is given by 



P 



[H,p]+C4p]. 



(4) 



The last term in Eq. Q describes incoherent processes 
such as spontaneous emission and is determined by 



'C7[p]=-7i(|3)(3|p-2|l)(l|p33 + p|3>(3|) 
-72(|3>(3|p-2|2)(2|p33 + p|3>(3|) 



(5) 



In a suitable interaction picture, the density matrix equa- 
tions follow as 



0-33 = - 2(71 + 72)0-33 + iG<J23 

- «G*0-32 - 1.9*0-31 , 
0-22 =272033 + iG*0-32 ~ «G'023 . 
1 



^5^13 



0-32 = 



0-31 = 



-(71 +72) +iA2 
iGo-22 - JGo-33 , 

^(71 +72) +iAi 



z.go-ii - 15033 , 



0-32 + tgai2 



0-31 + iGo-21 



(6a) 
(6b) 



(6c) 



(6d) 



0-21 = - [r - i(A2 - Ai)]o2i + iG*CT3i - zgo-23 ■ (6e) 

The detunings of the probe and the control fields from 
the respective transition frequencies are defined as 



Ai 
A2 



W31 - UJl , 
Ci;32 — U!2 ■ 



(7a) 
(7b) 



We have further included pure dephasing of the ground 
state coherence, e.g., due to phase changing collisions, 
and denote the total decay rate of the coherence by F. 
The remaining density matrix equations follow from the 
constraints an + 022 + 033 = 1 and o-y = a*^. 



B. Propagation equation 

We use Maxwell's wave equations to simulate the spa- 
tial evolution of the control and the probe beams through 
the medium, in order to study the effect of both diffrac- 
tion and dispersion during the propagation. The wave 
equations for the probe {j ~ 1) and control {j ~ 2) fields 
can be written as 



expressed in terms of both the atomic coherences as well 
as the susceptibility as 



(9) 



Vj = Af [dsjas.e-'^^' + c.c. 



[XsjBjt'je 



+ C.C.) , 



where Af is the density of the atomic medium, and X3i 
and X32 are the susceptibilities for the response to the 
probe and control fields, respectively. In slowly varying 
envelope and paraxial wave approximation, Eqs. ([IJ, ^ 
and ([HI) lead to propagation equations for the two fields 
given by 



dg I / 9^ 9^ , 
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(10a) 
(10b) 



The first terms in the parentheses on the right hand sides 
account for the diffraction. The second terms on the right 
hand sides are responsible for the dispersion and absorp- 
tion of the both the control and probe beams. Note that 
the two propagation equations are coupled via the sus- 
ceptibilities X31 and X32- 



C. Medium susceptibility 

Next we calculate the response of the medium to the 
probe and control fields, characterized by the respective 
susceptibilities. In steady state, the atomic coherences 
03i(wp) and 0-32 (cjc) are obtained from Eqs. ([6|) as 



N- 



03j = 



3j 

D 



(11) 



where the numerators iV32, N^i and the denominator D 
arc listed in Appendix|21 The expressions are rather com- 
plex, since we include the fields to all orders, in order to 
account for nonlinear effects. To simplify the expressions, 
we have assumed equal decay rates on the two transitions, 
71 = 72 = 7- Using Eqs. Q, one can readily obtain the 
susceptibility at the frequencies at Wc and ojp as 



, . A^M32P 
X32(Wc) = ^^ 0-32 , 

X3l(Wp) = — cr3i. 

D. Beam profiles 



(12a) 
(12b) 
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where Vj are the macroscopic polarizations induced by 
the control and probe fields, respectively. They can be 



In the main part of our result section, we will numer- 
ically propagate complex transverse beam profiles. But 
first, in order to interpret the effect of the beam profiles 
on the propagation, we chose the transverse spatial pro- 
file of the control field as a Hermite-Gaussian. At z = 0, 



it can be written as, 



Gix,y)^GoH„ 



V2\ ,, fyV2^ 



Wr 



H, 



(13) 



Here, Go is the input amplitude, and Wc is the width of 
the control field. The function Hk is a Hermite polyno- 
mial of order k, and the indices m and n determine the 
shape of the control field profile along the x and y direc- 
tions, respectively. Since wc want to consider the transfer 
of arbitrary spatial information, wc will study different 
values of TO, ri in the following. 

Similarly, the probe field is initially assumed to have a 
super- Gaussian transverse profile given by 



lix,y) =50 e 



— zra — 



(14) 



that this in turn leads to a smaller feature size of the 
image cloning scheme. 

The transparency of the medium to the probe field can 
be understood as arising from coherent population trap- 
ping (CPT) or electromagnetically induced transparency 
(EIT), depending on the relative strength of the two ap- 
plied fields. In both cases, destructive interference occurs 
between two different excitation pathways from |1) or |2) 
to the excited state |3). CPT and EIT are restricted to a 
certain transparency window. This transparency window 
can be controlled via the external parameters, and per- 
mits to transfer the transverse distribution of the control 
field onto the transmission profile of the probe field. This 
is the key mechanism of cloning the control field profile 
to the probe field. 

The corresponding real part of the susceptibility is 



The initial peak amplitude and the width of the probe 
field are denoted by go and Wp, respectively. Instead of 
choosing super Gaussian as an initial profile, the shape 
of the probe field can be consider any arbitrary shape 
such as a plane wave, Gaussian or hyperbolic shape. The 
desired spatial profile of the probe beam can be generated 
by using a spatial light modulator based on either liquid 
crystal or coherent EIT media[3lLl33|. 



III. RESULTS 
A. Linear and nonlinear susceptibility 

We now turn to our results, and start by studying the 
probe and control field susceptibilities at frequencies wi 
and UJ2 using Eqs. (|12[). For this purpose, we first consider 
a super Gaussian probe field and a Gaussian control field 
[m = 0, n = in Eq. ((T^ ] such that the probe suscepti- 
bility becomes inhomogeneous along the transverse direc- 
tions. Results are presented for two different cases of the 
initial field amplitudes. We denote the case go <C Go as 
the weak field limit, and go ^ Go as the strong field case. 
In our numerical calculations, we choose go = 0.015 7 and 
go = 0.157 for the two cases, respectively. 

Fig. [5] displays the spatial dependence of the probe 
field susceptibility. The different curves show Im(xi3) 
and Re(xi3) for three different detunings and both probe 
intensity cases as a function of x, with y = 0.005 cm. 
From panel (a), we find that the probe field is essen- 
tially transmitted without absorption in regions in which 
the control field amplitude is large, whereas in the wings 
(|a;| > 0.05 cm) with low control field, strong absorp- 
tion occurs. As expected; wc found that the width of 
the transparency window depends on the intensities of 
the control and probe fields as well as the decay rate of 
ground state atomic coherence. Comparing the two in- 
tensity cases, we find that a relative increase of the probe 
field intensity results in steeper transitions between re- 
gions of high and low absorption. Later, we will show 
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FIG. 2. (Color online) Probe beam susceptibility as a func- 
tion of the transverse coordinate x. The control field has a 
Gaussian (rn = 0, n = 0) beam profile, (a) shows the imag- 
inary and (b) the real component of the susceptibility. The 
three different curves show probe field detunings Ai = —0.005 
(blue dotted line, blue dashed line), Ai = 0.0 (green solid line, 
green dot-dashed line), and Ai = O.OO57 (red dashed line, red 
dashed double-dot line) for two data sets of probe amplitudes 
g = 0.0157 and g = O.I57, respectively. The other parame- 
ters are y = 0.005cm, uip = 150^m, lOc = 400/im, A2 = 0, 
r = O.OOI7, and N = 5xl0" atoms/cm^. 



shown in Fig. [HJb) . We notice that it can be controUcd 
between negative and positive values via the detuning. 
For red detuning of the probe field (Ai = 0.0057), it has 
a local maximum around x = 0, whereas for blue de- 
tuning, a minimum occurs around a; = 0. For non-zero 
detunings, the real part becomes spatially dependent in 
particular around the regions in which the absorption 
undergoes a transition from low to high absorption. 

Therefore, at red detuning, the transverse profile of the 
control field allows one to imprint a fiber-like refractive 
index gradient onto the atomic medium. This parabolic 
refractive index variation causes focusing of the probe 
field towards the center of the control field and also guides 
the probe field propagation along the propagation axis. 
Conversely, at blue detuning of the probe field, the real 
part has maxima at its wings, which is referred to as anti- 
waveguide like refractive index [33 . This anti- waveguide 
refractive index leads to shape distortions of the probe 
field. Hence, electromagnetically induced focusing and 
defocusing is possible in our setup by properly detuning 
the probe field. At two photon resonance condition, i.e., 
Ai = A2 = 0, only the imaginary part of the susceptibil- 
ity varies spatially, while the real part is constant. Then 
the probe field propagates through the transparency win- 
dow without being focused or defocused. Note that the 
probe field will nevertheless spread due to the inevitable 
diffraction. 

Next, we consider the effect on the control field. Fig- 
ure[3]shows the real part of the susceptibility experienced 
by the control field. It in particular illustrates that it be- 
comes especially important when the relative intensity of 
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FIG. 3. (Color online) Real part of the control field suscep- 
tibility as a function of transverse position x. The control 
field has a Gaussian (m = 0, n = 0) shape. The three dif- 
ferent curves show probe field detunings Ai — —0.005 (blue 
dotted line, blue dashed line), Ai = 0.0 (green solid line, 
green dot-dashed line), and Ai = 0.0057 (red dashed line, red 
dashed double-dot line) for two data sets of probe amplitudes 
g = 0.0157 and g = O.I57, respectively, with different scales 
of the axes. The other parameters are same as in Fig((2]). 
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FIG. 4. (Color online) Probe field susceptibility as a function 
of the transverse spatial coordinate x. The control field has 
a doubly-peaked spatial profile characterized by m = 1 and 
n = in Eq. (|13|l . The central part around a; = is magnified 
in the inset in the top part of the figure. The other parameters 
are y — 0.005 cm, Go — 1-57, go = O.OI57, Wp = 400/im and 
N = 7.5x10" atoms/cm^. 



the probe and control fields are of comparable strength, 
as for the intense probe field case, the spatial dispersion of 
the control field is about three orders of magnitude larger 
than in the weak probe field case. Hence, it is important 
to include atomic coherence effects in the paraxial wave 
equation for the control field. It is also evident from 
Fig. [3] that in addition to the diffraction, the variation 
of the refractive index induces further distortions of the 
control field shape for both red- and blue-detuned probe 
fields. Our further analysis also showed that the absorp- 
tion of the control field is rather low as compared to the 
probe field absorption in both cases, such that the con- 
trol field can be propagated through the medium without 
significant absorption for relevant parameter ranges. 

Finally, we analyze the dependence of the probe sus- 
ceptibility on the control field transverse beam shape. 
Figure m shows the transverse variation of the probe field 
susceptibility for a doubly peaked transverse spatial dis- 
tribution of the control field obtained by setting m = 1 
and n = in Eq. (jl3|) . The central part of the figure 
is shown enlarged in the inset. Overall, the interpreta- 
tion is similar to the case of a singly-peaked control field 
beam profile. But the doubly-peaked structure leads to 
the formation of two transparency windows. For positive 
probe detuning (Ai = 0.005), each of the transparency 
windows features a fiber-like refractive index, focusing 
the probe field towards the waist of the control field and 
thereby reducing the feature size. The generalization of 




FIG. 5. (Color online) Normalized intensity profile of the 
propagating probe beam against the transverse coordinate x 
for J/ = 0. (a) shows the beam profile at different propagation 
distances z. In (b), the transmitted probe beam intensity 
at the output of a 4 cm long medium is shown for different 
detunings. All other parameters are chosen as in Fig. [2] 



approximately z ~ 2 cm. Interestingly, the probe field 
width at this propagation distance is about 60/xm which 
is an order of magnitude less than the initial width of the 
control beam. At this distance, the Rayleigh length of 
the cloned probe beam is about 1.42 cm. Consequently, 
the cloned beam suffers from diffraction throughout the 
further propagation. We also notice from Fig. [5] that 
at two photon resonance condition Ai = A2 = 0, the 
probe beam experiences diffraction. In contrast, at red 
detuning, the probe field propagates unaltered through 
the medium. We thus find that the probe field diffraction 
can be controlled by a proper choice of the single photon 
detuning Ai of the probe field. 

Fig. ini^b) further illustrates how the intensity and de- 
tuning of the probe field can be used to control the width 
of the transmitted probe beam. It shows that the probe 
beam is more tightly focused by a factor of about two for 
the strong control field case compared to the weak probe 
field case. At red detuning condition, our numerical com- 
putations show that the width and the transmission of 
the output probe beam are 85/^m [55/^m] and 20% [5%] 
for weak [strong] field case. 

Figure [5] show the corresponding results for the spa- 
tial variation of the transmitted control beam intensity 
profile in the y = plane after propagation through a 4 
cm long atomic medium. For the weak probe field case, 
the control field propagates essentially without any dis- 
tortion. But in the strong probe field case, we find from 
Fig.[6]that the control field shape is distorted. Neverthe- 
less, the integrated transmission intensity of the control 
field is approximately 100% in both cases. Thus, we con- 
clude that the control beam intensity is spatially redis- 
tributed in the strong probe field case. As a consequence, 
it is necessary to include the effect of atomic coherence in 
the propagation of the control field in the non-resonant 
as well as in the resonant condition for the strong probe 



these results to other spatial modes of the control field 
with different values of m, n is straightforward. 



B. Propagation dynamics of probe and control 
beams 



In this Section, we present results for the evolution 
of the spatial beam profile throughout the propagation. 
For this, we numerically integrate the full set of parax- 
ial wave equations Eqs. (jlOp using a higher order split 
operator method [30|. We begin with a super-Gaussian 
probe beam shape and a Gaussian control beam shape to 
demonstrate the mapping of the control field profile onto 
the probe field. Results are shown in Fig. [SJ Panel (a) 
depicts the intensity profile of the probe beam at differ- 
ent propagation distances z. The initial field amplitudes 
are Go = I7 and qq = O.I57, corresponding to the strong 
probe field case. We find that the probe beam profile 
gradually acquires the control beam shape as it propa- 
gates through the medium. This process is completed at 
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FIG. 6. (Color online) Spatial intensity profile of the con- 
trol field as a function of the transverse coordinate x after 
propagation through a 4 cm long atomic medium. Results 
are shown for different field parameters. The profile is shown 
in the y = Q plane. All other parameters are the same as in 
Fig. El 
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FIG. 7. (Color online) Transverse spatial intensity profile 
of the probe and control fields after propagation through an 
atomic vapor cell of length 4cni. The initial control field pro- 
file is chosen as a Hermite- Gaussian doubly-peaked profile. 
The top panel (a) shows the probe field, the bottom panel 
(b) for the control field. The parameters are as in Fig. |4] 
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FIG. 8. Transverse spatial intensity profile of the control and 
probe fields, (a) shows the initial profile of the control field 
at the medium entry, (b) shows the control field profile after 
propagation through a vapor cell of length 3cm filled with 
rubidium atoms, (c) shows the probe field profile after the 
propagation. Initially, the probe field is chosen as a plane 
wave. The parameters are g = 0.157, G = 1.67, A2 — 0, 
Ai = 0.0027, r = O.OOI7, and iV = 10^^ atoms/cm^*. 



field case. 

Next, we study cloning of a two-pcakcd Hcrniitc- 
Gaussian control beam onto the spatial profile of the 
probe beam. Fig. [7] shows the peak-normalized inten- 
sity of both the transmitted probe and control beams 
after a propagation distance of 4 cm. It can be seen 
from Fig. [7] that while the overall structure of the control 



beam is preserved, the width of the cloned probe beam 
profile is reduced by a factor of about 2.5 as compared 
to the initial width of the control beam. Therefore, the 
finesse of the transmitted probe beam is about 2.5 times 
greater than that of the control beam, which is consis- 
tent with the findings of an recent experiment by Li et al 
[23| . The integrated transmission intensity of the cloned 
output probe beam is 5% of the integrated intensity of 
the input probe beam. The transmission of the cloned 
beam can be increased by decreasing the optical density. 
We also found that the width of the cloned beam can 
be decreased by decreasing the relative intensity of the 
probe and control field giving us control over the finesse 
of the cloned probe beam. 

Finally, Fig. |S] shows how a complex image encoded in 
the transverse control field intensity profile can efficiently 
be cloned onto the probe field. For this, we choose the 
three letters "CPT" as initial profile of the control beam, 
and a plane wave for the initial probe beam profile. We 
find from Fig.[Sl[b) that the control beam profile is signif- 
icantly distorted after 3 cm propagation length through 
a vapor cell. This distortion arises from diffraction and 
from the variation of the refractive index experienced by 
the control field. In contrast, the fiber-like spatial dis- 
persion of the probe field can be used to compensate or 
to reverse the effect of the diffraction. As a result, the 
cloned image in the transverse profile of the probe beam 
is tightly focused. We further note from panel (c) that 
the feature size in the cloned image is reduced by a fac- 
tor of 2 as compared to the incident control field image. 
This suggests that the finesse of the transmitted cloned 
image is two times greater then that of the incident con- 
trol image. Hence, the feature size of the cloned image is 
decreased twofold. 



IV. CONCLUSION 

In conclusion, we have studied the possibility of cloning 
of an arbitrary images encoded in the spatial profile of 
a control beam onto a probe beam. Our method is 
based on coherent population trapping in a three-level 
atomic lambda system. We have considered both, weak 
and strong probe fields with strength comparable to the 
control field, and have discussed the differences for light 
propagation in the weak and strong field limits. In the 
"strong field limit" where the probe and control fields are 
of comparable strength, we firstly calculated the atomic 
susceptibilities including both linear and nonlinear effects 
for the two fields. We then found that a waveguide-like 
structure can be formed inside the medium at red de- 
tuning of the probe field when applying a Gaussian con- 
trol and a super-Gaussian probe. At the same time, a 
transparency window centered in the waveguide, which 
can be controlled by changing the relative intensities of 
the probe and control fields. This transparency window 
allows to transfer the transverse intensity profile of the 
control onto the transmitted probe field. By numeri- 



cally solving the propagation equations for both fields, 
we found that the spatial profile of the control is grad- 
ually mapped onto the transmitted probe already after 
a few Rayleigh lengths. Interestingly, the feature size in 
the spatial profile of the probe field is reduced compared 
to that of the original control field structure. In order 
to show that our method works for arbitrary images, dif- 
ferent spatial profiles of the control and probe fields are 
considered. In particular, the three letters "CPT" ini- 
tially encoded on the control field are cloned onto the 
transmission profile of the probe. In this process, the fea- 
ture size is decreased by a factor of 2. Interestingly, even 
though cloning of arbitrary images onto the probe field 
is constructed nicely, the control field image is severely 

I 



distorted due to diffraction throughout the propagation. 
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Appendix A: Coefficients for susceptibility 



^31 =(|Gp(7(*7 + Ai)(r2 + (A2 - Ai)2) + (7(ir + A2 - Ai) 
^32 =(|ff|'(7(«7 + A2)(r2 + (A2 - Ai)2) + 7(tr - A2 + Ai)|.g, 
D =7|G|S + IGI"* [3|5p(7 + 2r) + 27(7r + Ai(A2 - Ai))] 

+ 7|<?P [2|5p(7r + A2(Ai - A2)) + \g\' + (7' + A2) ((A2 
+ |Gni.gp ((47 + r) A^ + 27r(27 + BE) + 2(r - 47) A2 Ai 
+ 3|5|^(7 + 2F) + 7(7' + A?)((A2 - Ai)^ + F^)] . 



+ F(A2 + Ai))|.g|2 + 7(zr + A2 - A,)\G\^)g) (AI) 
' + (7(ir-A2Ai)+r(A2 



A,))\Gf)G) (A2) 



Ai)2+f2)] 
(47 + F)A?) 



(A3) 
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